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Abstract 

We show the equivalence between two different communication schemes 
that employ a couple of modes of the electromagnetic field. One scheme 
uses unconventional heterodyne detection, with correlated signal and im- 
age band modes in a twin-beam state from parametric downconversion. 
The other scheme is realized through a complex-number coding over quadr- 
ature-squeezed states of two uncorrelated modes, each detected by ordi- 
nary homodyning. This equivalence concerns all the stages of the commu- 
nication channel: the encoded state, the optimal amplifier for the chan- 
nel, the master equation modeling the loss, and the output measurement 
scheme. The unitary transformation that connects the two communica- 
tion schemes is realized by a frequency conversion device. 

1 Introduction 

The fundamental theorems of quantum communication theory ^] establish 
an upper bound on the mutual information for all possible coded inputs, and 
all kinds of output detection at the channel. These results, when applied to 
a single narrow-band linear bosonic channel ||, specify the ultimate quantum 
capacity per use (i.e. the maximum mutual information per mode subjected 
to any constraint on the channel), which, for fixed average-power constraint, 
could be ideally achieved by direct detection of number-states with thermal a 
priori input probability. Any other communication scheme can transmit just 
less information than the number-state (NS) channel. At high power levels, for 
example, the ideal NS channel achieves 1.44 bits per use more than the ideal 
coherent-state (CS) channel, and 0.44 bits more than the quadrature-squeezed- 
state (QS) channel ||. (The CS channel utilizes heterodyne detection of CS, 
with Gaussian a priori probability; the QS channel uses homodyne detection of 
QS having Gaussian probability and optimal fraction of squeezing photons). 
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Although the NS channel is the optimal one, several reasons relegate it to 
a purely theoretical plan, and lead to consider the QS channel as a more re- 
alistic (and not much less efficient) communication scheme. First of all, it is 
a challenge to achieve number eigenstates experimentally. Second, the optimal 
amplification for the QS communication channel is achieved by a phase-sensitive 
ideal amplifier, which can be realized experimentally, whereas a concrete real- 
ization of the ideal photon number amplifier for the NS channel is still unknown 
[1,0- Finally, the detrimental effect of loss, widely neglected in the literature 
[5 1, easily degrades the NS channel: for example, for ten photons of average 
power, a signal attenuation of 0.2dB is enough to degrade the NS channel under 
the capacity of the CS one, and for higher power the effect is more dramatic 
@. The QS channel is more robust to losses than the NS channel, and at low 
powers and not too high losses it remains above the capacity level of the CS 
channel 0. 

In this paper we show the equivalence between two communication schemes 
that employ complex-number encoding carried by two-mode states. The poten- 
tiality of a quantum communication channel that uses two modes has already 
been noticed, in particular for phase-modulation-based digital communications 
The two equivalent schemes are the following: 

(i) a complex number is encoded on a twin-beam state generated by parametric 

downconversion; decoding is achieved through unconventional heterodyne 
detection of both the signal and image-band modes which form the corre- 
lated twin-beam; 

(ii) the real and the imaginary part of the complex number are independently 

encoded over two uncorrelated QS pertaining two different modes; the 
two QS are decoded through ordinary homodyne detection on each mode 
separately. 

The equivalence between the above schemes involves all the stages of the 
communication channel: coding, decoding and optimal amplification. More- 
over, it includes the master equation that models the loss along the line. The 
equivalence is realized by a unitary transformation that physically corresponds 
to 50-50 frequency conversion. The equivalence between the two channels is 
relevant for applications, because the heterodyne twin-beam communication 
channel is easier to achieve experimentally than the QS channel. 

The paper is organized as follows. In Sec. II we describe the twin-beam 
communication scheme. We introduce the eigenstates of the heterodyne de- 
tector photocurrent, and give a scheme that realizes states that approach such 
eigenstates used to encode the transmitted information. In Sec. Ill we discuss 
the two-mode master equation that models the effect of both losses and linear 
phase-insensitive amplification, and then we derive a Fokker-Planck equation 
for the unconventional-heterodyne probability. In Sec. IV the channel equiva- 
lence is analyzed. We show that the 50-50 frequency conversion disentangles 
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a twin-beam state into a couple of QS, whereas, at the same time, it changes 
a phase-insensitive amplifier (PIA) into two independent phase-sensitive ampli- 
fiers (PSA). On the other hand, the master equation for the loss turns out to 
be invariant under frequency conversion. In Sec. V the Fokker-Planck equations 
for the two schemes are compared, and the mutual information that character- 
izes the communication channels is evaluated. Finally, maximization over input 
probabilities is carried out to get the capacity per use. 



2 The twin-beam communication scheme 
2.1 Unconventional-heterodyne detection 

Under ideal conditions, the output photocurrent of a heterodyne detector cor- 
responds |l0|, |n] to the complex operator Z = a + b\ a and b denoting the 
annihilator of the signal and the image-band mode, respectively. Ordinary het- 
erodyne detection corresponds to the joint measurement of both quadratures 
of the field a, with b as the vacuum image-band mode responsible of the 3dB 
heterodyne added noise. The communication scheme here presented is based on 
unconventional heterodyne detection, namely with the signal and image-band 
modes both nonvacuum. The eigenvector of Z with complex eigenvalue z can 
be written in several equivalent forms. In terms of the eigenvectors \x)<j, of the 
quadrature = ^(c l e 1 ^ + h.c.) of the mode c — a,b one has JO, |l2] ] 

^e 2lxlmz \x) ® |Rez-x> 

V 7T 



— OO 

+00 



^Le- 2i v Rez \y + Imz) w/2 ®\y) n/2 , (1) 

where \ip) £g> \(p) denotes a vector in the two- mode Hilbert space H — 7i a ® Tib 
[the notation | )) remembers that the state is a two- mode one]. A number 
representation of the state ([[]) can be found in Ref . |Q . The eigenstate for zero 
eigenvalue reads 

1 1 00 

|0» = -= exp(-at&t)|0) ® |0> = -= ]TR» ® |n) . (2) 

It is convenient to write the eigenstate of Z corresponding to complex eigenvalue 
z as follows 

\z)) = L^ ex p(- a W)\z)®\%, (3) 

V 7 !" 

where \z) denotes a customary coherent state. The states {\z})} are Dirac- 
normalized as ((z\z')) = 8^ 2 \z — z'), and form a complete orthogonal set for 
H. 
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In the following, we will introduce a set of physical (normalizable) two-mode 
states that approach the (infinite-energy) eigenstates of Z. 



2.2 Limited power near-eigenstates of the photocurrent 

The physical realization of states that approach the eigenstates \z)) in Eq. (||) 
is suggested by recognizing that the zero-eigenvalue state (Eh is just a "twin- 
beam" at the output of a PIA in the limit of infinite gain . All the other 
states with z ^ can be approached by suitably displacing cither a single mode 
a or b, or both of them, and this can be done by means of high-transmissivity 
beam-splitters with a strong coherent local oscillator. Such physical realizations 
\z))\ can be written as follows 

\z)) x = D a (v)D b (w)e t!inh ~ lx< - ab ~ a>b ^\0)®\0) (z = v + w), (4) 

where D c (u) = e uct_uc denotes the displacement operator of the mode c = a,b, 
and e tanh lA ( ab - atfct ) j s the unitary transformation achieved by a PIA with 
gain G = (1 — A 2 ) -1 . [The PIA Hamiltonian (a)b^ + h.c.) corresponds to a 
parametric downconversion in the rotating- wave approximation and for classical 
undepleted pump] . The average number N of photons per mode of the state ([|) 
is given by 

The probability density of the photocurrent in the state \u)}\ has the Gaussian 
form 

1 



l«^»A| 2 = rXTexp(-^-r^) , (6) 



with variance 



A 2 = (l-A)/(1 + A). (7) 

Using the identity 

e -«ai,-at6t) D a (v)D b (w)e^ ab - aib ^ = D a (vcbZ + wshZ)D b {vsbt + wchO (8) 
Eq. (0) rewrites as follows 



gtanh 1 X (ab— ) 



Xw 



VT^A 2 



Xv 



VT^A 1 



(z = v + w) . (9) 



Hence, the state \z}}\ can be also obtained through phase- insensitive amplifi- 
cation of input signal and idler which have the precise amplitude relation in 
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Eq. (|j). Notice that the probability density (||) does not depend explicitly on 
v and w (the share of displacement of the two modes), while the mean number 
of photons does. The constraint v + w — z over states (J4j) and (^) implies that 
for each z there is a family of states approaching the eigenstate \z)) according 
to Eq. (j6|). For v — w — z/2, the most symmetrical state 

\ Z )) X = gtanh- 1 A(a 6 -at h t) g \ 2 /2A X ) (10) 

achieves the best phase sensitivity, as shown in Ref. [[12]. For high photon 
number, the marginal probability density for the phase <f> = arg(Z) is Gaussian, 
and in the limit of infinite PIA gain (A — > l - ) the r.m.s. phase sensitivity is 
optimized versus the total photon number n = 2N to the value H] 

= (A<j> 2 )^ ~ * • (11) 
V2 n 

For a repeatable phase measurement scheme on a two-mode field, see Ref. jD| . 

The couple of coherent states on the right of Eq. ( |l0| ) can be generated by 
means of a single coherent state and a frequency conversion device with suitable 
pump strength and phase. Indeed, for any complex a, one has 

exp [J (e iargQ a+& - e _<a ' 80t a6 + ) |V5a) O |0) = |a) ® |a) , (12) 

and the unitary operator on the left side of Eq. ( |l2[) describes a frequency con- 
version device [in the parametric approximation of classical undepleted pump, 
this can be realized through a three-wave (or degenerate four-wave) mixing in 
a nonlinear (x' 3 ^) medium]. The two equivalent experimental set-ups to 
generate the state ( |l(i| ) are sketched in Fig. 1. 



3 Losses and distributed amplification 

The effect of losses on the communication channel can be modeled by the fol- 
lowing master equation 

d t Q = C r g = T{(n a + l)L[a] + (n b + l)L[b] + n a L[a)] + n b L[tf]} g . (13) 

In Eq. (|l3|), the superoperator Cr gives the time derivative of the density 
matrix g of the radiation state in the interaction picture, and acts on g through 
the Lindblad superoperators L[c]g — cgc^ — ^(c^cg + gc^c) Jig ]. The damping 
rate F is supposed to be equal for both modes, whereas the mean number of 
thermal photons n a and rib at the frequency of modes a and b can be neglected 
at optical frequencies. The absence of cross-terms that correlate the two modes 
is a consequence of the rotating-wave approximation assumed in the ordinary 
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derivation of the master equation Jl7[ . In a similar fashion, an active medium 
amplifier in the linear regime can be described by the superoperator 

C A =A{(m a + l)L[J] + (m b + l)L[b t }+m a L[a]+m b L[b]} , (14) 

where A denotes the gain per unit time — i. e. the amplifier length — and m a 
and m b are related to the population inversion of the lasing levels at resonance 
with a and b, respectively. Finally, a distributed parametric amplification with 
modes a and b correlated by a classical pump is represented by the following 
commutator 

C K = K[a^tf - ab, •] , (15) 

where K is the gain per unit time and is related to the intensity of the pump. 
The following differential representation of modes 

4*» = (f-0*)l*» a t l4 = (| + ^)l*» 

b\z)) = (^-d z )\z)) 6 f |4 = (f + (16) 

converts a two-mode master equation of the general form 

Cg t = 2{(A + C a )L[a^ + (A + C b )L[b^} + (B + C a )L[a} + (B + C b )L[b}}g t 
+ K[a^tf-ab,Q t } (17) 

into the following Fokker-Planck equation 

d t P{z,z;t) = {Q {d z z + d,z) + 2D0 2 Z ,} P(z,z;t) , (18) 

where P(z, z; t) = ((z\g t \z)) denotes the (unconventional) heterodyne probability 
density, and the drift and diffusion terms are given by Q = B — A — K and 
D = A + B + C a + C b . Notice that the coefficients in the master equation 
( p~7[ ) are not independent (the difference of the first two coefficients equals the 
difference of the last two in the curly brackets), and this is the condition under 
which all the derivatives in Eq. ([l6]) gather to give the simple Fokker-Planck 
equation (H). Hence, when considering the effect of loss in Eq. (|l3|), the 
assumption of equal damping for the two modes is crucial. Of course, the same 



argument holds true for the gain A in Eq. (14). 

The solution of Eq. ( |T^ ) will be given in Sec. V, where also the mutual 
information transmitted by the channel will be evaluated. In the next section 
we will show the equivalence of the twin-beam channel with a communication 
scheme based on a couple of uncorrelated quadrature-squeezed states. 
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4 Equivalence between squeezed-state and twin- 
beam channels 



Let us consider the unitary transformation that describes a 50-50 frequency 
conversion device from mode a to mode b. Under suitable choice of phases, the 
corresponding unitary operator U writes 



U = exp 



(19) 



so that the Heisenberg evolution of modes is given in matrix form as follows 



u 



1 

V2 



1 -1 
1 1 



(20) 



The action of the operator (|T^) on a twin-beam state is more easily evaluated 
using Eq. ([I|) for the twin-beam with v = w = z/2. Since the vacuum state |0) <8> 
|0) is an eigenvector of the frequency conversion operator ( |l9| ) with eigenvalue 
1, from the identity 



e A Xe- A =X + [A,X} + -{A, [A,X]} 



one has 

U\z))x 



zlmz \ 

Rez)_ A 



5 a (tanh" 1 A) 



|0) ® 



/Rez 
W2 



j 5 b (-tanh _1 A) 



D a 

= \ilmz)x 

where S c (r) is the squeezing operator for the mode c = a, 6, namely 



— i c' — c 2 



r 

S c (r) = exp y- 

and the squeezed state \a)\ is defined as follows 

\a) x = D c (a/V2) ^(tanh^ 1 A)|0) 



(21) 

10) 
(22) 

(23) 

(24) 



Hence, by means of frequency conversion the twin-beam state (||) disentangles 
into two squeezed states which are still related in intensity and phase. [For 
modes a and b at the same frequency and different wave vectors or polarization, 
this "disentanglement" can be also achieved by means of a 50-50 beam splitter 
[[L8), i.e. by a passive device]. 

As regards the effect of frequency conversion at the output of a lossy / amplified 
channel, notice that a superoperator of the form 



C = aL[a] + pL[b] + 7 L[a t ] + 6L[tf] + C 



K 



(25) 
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undergoes the following transformation 

UCfr = ^±l(L[a} + L[b}) + 2±^(L^}+L[b^)+^[a 2 -^\-} 

K 2 + 2 . a — j3 . . 7 — (5 . . 
— — [o — o 1 , -J H — (cross-terms) H — (cross-terms) . (26) 



We are not interested in the cross-terms in Eq. ( p6[ ) that act jointly on H a <&7i.b, 
because they vanish when n a = n b and m a — m& in Eqs. (|T|) and @, 
respectively. In this case one has 

U£ r U 1 = C r ; UC A U^ = £ A . (27) 

Hence, the 50-50 conversion leaves the superoperators Cr and £a (for both loss 
and PI A) invariant. This means that the disentanglement of the twin-beam 
occurs equivalently at whatever time during transmission. On the contrary, 
a distributed PIA with pump-correlated a and b is not invariant under the 
transformation (^0|). Indeed, a distributed PIA followed by frequency conversion 
is equivalent to a couple of independent PSA's, as shown by the commutator 
terms in Eq. (Eq), namely 



UCkU 1 = C'k = f [o 2 - « f2 , •] - f [b 2 6 f2 , •] (28) 

Notice that for K > the amplified quadrature components of the fields a and 
b through these commutators are Y a = and Xf, = 6q, respectively, which 
are the right quadratures in order to enhance the signal carried by the couple 
of squeezed states (p2|). Of course, the communication scheme that encodes the 
information on states ( |22|) needs two independent homodyne measurements of 
the quadratures Y a and Xf,. 

We have shown the equivalence of the two channels as regards the input 
states and the evolution master equation. It remains to show that also the final 
detection stage is equivalent in the two schemes. The unconventional heterodyne 
detection is described by the following orthogonal resolution of the identity 

dp,(z, z) = d 2 z S (2) {Z-z) = d 2 z \z)) ((z\ . (29) 



The unitary operator (19) transforms the orthogonal resolution fl29| ) as follows 

Udfl(z, z)^ = d 2 z S(V2X b - Rcz) 5(V2Y a - Imz) 

= dx dy 5(X b - x) 6(Y a - y) , z = V2(x + iy) . (30) 



The last orthogonal resolution in Eq. (|30j) is just the one corresponding to two 
independent homodyne measurements of quadratures X b and Y a . 

In conclusion of this section, we also show the equivalence between the QS 
scheme and the twin-beam scheme at the level of Fokker-Planck equations. Us- 
ing the homodyne probability density 

P a (y;t)P b (x;t) =Tr[&|y) f f (y| O \x) 00 {x\] (31) 
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the master equation 

C'et = 2 {A(L[at] + L[tf]) + B(L[a] + L[b})} g t 

+ |[ a 2- a t 2 ,^]_|[ 6 2_ 6 t 2 ^ t] (32) 

can be written in the Fokker-Planck form 

d t {P a (y; t)P b (x; t)} = |q (d x x + d v y) + j (d 2 xx + d 2 yy ) j P a (y; t)P b (x; t) (33) 

with Q = B - A - K and D = A + B. In obtaining Eq. (||) we used the Wigner 
rappresentation of both a and b modes, and then evaluated the marginal inte- 
gration over the quadratures X a and Y b . The equivalence of the Fokker-Planck 
equation (133) with Eq. (p~8) is evident, after the coordinate transformation 



z = V2{x +iy), z = V2(x - iy) , (34) 

and upon renaming the product of probabilities as follows 



P'(z,z-t) = -P a [-=;t )P b [-=■ 



t) . (35) 



The equivalence between the two communication channels is schematized in 
Fig. 2. 



5 Evaluation of the channel capacity 



The solution of the Fokker-Planck equation ( |l8|) for a Gaussian initial probabil- 
ity density 



P(z,z;0) 



1 



ttA 2 (0) 



exp 



\z — w\' 
A 2 (0) 



keeps the Gaussian form at all the time, and writes 



P(z,z;t) = 



1 



7rA 2 (i) 



exp 



we 



-Qt\2 



A 2 W 



where the "evolved" variance at time t is given by 

A 2 (t) = ^(1 - e~ 2Qt ) + A 2 (0)e- 2Qt 



(36) 



(37) 



(38) 



For the twin-beam communication scheme one has A 2 (0) = A|, according to 
Eq. (||). For the scheme based on homodyne detection over the couple of QS 
states (E2), one equivalently finds the Gaussian density 



o(x\®*(y\U\w)) x 



exp 



Row V 



Imw V 

717 



(39) 
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The mutual information of a communication channel with conditional prob- 
ability density p(w\z), encoded complex variable w, and a priori distribution 
p(w) writes as follows 

1= [d 2 w [ d 2 zp(w)p(w\z) In . (40) 

J J J a z w' p(w')p(w'\z) 

For a priori and conditional probability densities both Gaussian, i.e. 

p(w) = J^expf-MT) (41) 



one has 



Z = ln(l + ^). (43) 

According to the Eq. (^) for the average number N of photons per mode, the 
variance of the prior distribution for the twin-beam scheme under the constraint 
of fixed average power reads |l9| 

^ = 4^-^). (44) 

Hence, the mutual information transmitted at time t for a twin-beam channel 
modeled by the Fokker-Planck equation (|l^) is given by 



(45) 



Eq. ( [45| ) is valid also for the QS scheme, with parameter A and variance A| 
related to the squeezing parameter through Eqs. ( |22| ) and ([?]). In the lossless 
case, the mutual information is maximized by A = Nj (N + 1) to the value 

7 = 21n(l + 2iV) . (46) 

This information is always greater than the information of two lossless coherent- 
state channels, with one additional bit per mode for high average power (N 3> 1). 
For a lossy channel with damping rate T and n a = ri\, = n thermal photons one 
has D = T{2n + l)/2 and Q = T/2. In this case, for A = N/(N + 1), Eq. @ 
rewrites 



4N (N + 1 
1 + (2N + 1) (2n + 1) (e r * - 1) 
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The exponential erasure of information in Eq. 
master equation (p7|), with K = T/2, namely 

AN (N + l) 



(p7|) becomes linear using the 



I = In 1 



l + (27V+l)(2n + l)rt 



(48) 



This case represents an ideal distributed parametric amplification that works 
against the detrimental effect of loss. Notice that, due to the energy of the 
pump, the condition of fixed average power is not strictly satisfied. Indeed, 
after a transient, the increase of the mean number of photons is approximately 
linear in time as follows 



N(t) = \ 



N(0) 



n — — 



N(0) 



5" 



2rt 



+ T(2n+ l)t 



(49) 



6 Conclusions 

The maximum mutual information per mode for narrow-band linear bosonic 
channels under the constraint of fixed average power is achieved by the ideal NS 
channel with thermal input probabilities. Nevertheless, the experimental real- 
ization of number eigenstates and of an ideal photon number amplifier that can 
assist the channel is still unknown. Hence, our interest is turned into QS chan- 
nels, which are feasible and still lead to a satisfactory efficiency. In this paper 
we have shown a new communication scheme that is based on unconventional- 
heterodyne detection on two-mode states ("twin-beam") and represents an al- 
ternative way to achieve the QS channel capacity. We have given two experi- 
mental set-ups to obtain the twin-beam states. We have shown that the twin- 
beam communication channel is equivalent to a channel that employs a couple 
of QS channels and ordinary homodyne detection as decoding. The correspon- 
dence between the two schemes holds also in the presence of loss and optimal 
amplification. The equivalence is realized by a unitary transformation that 
physically corresponds to the 50-50 frequency conversion between the two field 
modes. The twin-beam scheme is easier to obtain experimentally as compared 
to the QS scheme. Both the encoding and decoding stages are simpler for the 
twin-beam channel than for the QS one. As encoder one just needs parametric 
downconversion and coherent states, instead of squeezing. On the other hand, 
as decoder one has just one heterodyne detector versus a delicate balancing of 
two equal homodyne detectors. As compared to the single homodyne QS chan- 
nel, the heterodyne one, by employing two field modes, carries a double amount 
of information. 
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Figure 1: Outline of the two alternative experimental set-ups to generate 
the two- mode states \z))\ of Eq. (|Io|). The labels H and LO denote high- 
transmissivity beam splitters and local oscillators respectively, which realize the 
displacement operators in Eq. (^J) . The input states at the parametric amplifier 
(PIA) are the vacuum state in the upper scheme, and two independent coherent 
states in the lower one. In the last case the couple of coherent states is generated 
by a single coherent state through frequency conversion (FC). 
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Figure 2: Equivalence between the channels transmitting twin-beam states \z))\ 
and couples of quadrature-squeezed states |zlmz),\ ® |Rez)_A- The decoding 
measurements are the unconventional-heterodyne detection (complex photocur- 
rent I z ) and two independent ordinary homodyne detections (real photocurrents 
I x and I y ), respectively. The superoperators Cr, Ck and C' K model the loss in 
Eq. ( |l3| ) , the phase-insensitive amplification in Eq. ( |l5| ) and the phase-sensitive 
amplification in Eq. (|2^), respectively. The equivalence between the two com- 
munication schemes is realized by the unitary transformation (pof), namely a 
50-50 frequency conversion. 



15 



LO 



oo 
o> 
o\ 

c 

00 

(N 
> 

m 
o 
^- 
o 

a,: 

CT 1 - 



10) 
|0) 



PIA 

gain(l - A 2 )" 1 



"twin-beams" 

|0»A 



pump 




z/2A x ) 
z/2A x ) 



pump 



H 



H 



LO 



pump 



a 


PIA 




— ► 






gain(l - A 2 ) -1 









X 



Fig. 1. 

G. M. D'Ariano and M. F. Sacchi 

Equivalence between ... 



1 



oo 
o> 
o\ 

c 

00 

(N 
> 

m 
o 
^- 
o 

a,: 

CT 1 - 



\ilmz)\ 

|Rez)_ A — 



C r + C 



A 



A 



het. det. I 7 oc a + 



horn. det. I y oc (a — a j( )/2i 
horn. det. J x . oc (b + 6 f )/2 



X 



Fig. 2. 

G. M. D'Ariano and M. F. Sacchi 

Equivalence between ... 



1 



